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C^- , Abstract 

VO ■ 

£Sj ■ We examine the prospects for using lead as a supernova-neutrino detec- 

tor by considering the spectrum of electrons produced, and the number of 
one- and two-neutron events. We show that the electron energy spectrum 
from charged-current reactions can be used to extract information about the 
high-temperature component of the neutrino spectrum. Some degree of elec- 
i-£h I tron neutrino oscillation is expected in the supernova envelope. We examine 

the prospects for untangling the signatures of various oscillation scenarios, in- 
Q_i| eluding, e.g. normal or inverted hierarchies, and different values for the small 

mixing angle #13. 
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■ 1 Introduction 

The idea of detecting supernova neutrinos is exciting because although galactic supernovae 
are rare, we could potentially learn a great deal from the neutrinos. Measuring their 
spectra would give us information about the mass of the proto-neutron star core and its 
equation of state, and provide input for supernova explosion calculations. In addition, 
we can add to what we've learned about neutrino oscillations from solar, atmospheric, 
reactor, and accelerator neutrinos. In this paper, we analyze the capabilities of a detector 
based on lead Jl], || that can see both electrons from charged-current neutrino interactions 
and neutrons from charged- or neutral-current interactions, e.g. a detector containing lead 
perchlorate ||. 

Supernova-neutrino detection differs from solar-neutrino detection. Supernova neu- 
trinos have a higher average energy, and an intense flux for a very short period, ~ 10 s. 
Calculations of neutrino diffusion in the proto-neutron star show that neutrinos emitted 
as /i and r types have higher average energies, (E u ) = 16 — 25MeV, than those emit- 
ted as electron or anti-electron types, which have {E Ue ) ~ HMeV and (Ep e ) ~ 13MeV 
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respectively. The energy distributions can be fit to Fermi-Dirac spectra, although some 
calculations show that they differ from Fermi-Dirac form on their tails Q. This has been 
discussed extensively in ||. The interesting supernova physics conveyed by neutrinos lies 
in the details of the energy distributions, acquired as the neutrinos are emitted from the 
supernova core. Neutrino oscillations will mix the spectra in the outer envelope of the 
supernova. 

Recent data from the Sudbury Neutrino Observatory (SNO) and Superkamiokande 
(SuperK) have demonstrated that neutrinos oscillate, implying that they have mass. Su- 
perK indicates that atmospheric n neutrinos oscillate into objects that are not electron 
neutrinos, with <5m 2 tmos « 3 x 10~ 3 eV 2 and a mixing angle of sin 2 2# atm0 s ~ 1 |§- Al- 
though it has long been known that the flux of electron neutrinos from the sun is smaller 
than expected, SNO has very recently used its sensitivity to neutral-current scattering 
to measure the total flux from all species of neutrinos. The flux is approximately the 
same as predicted by the standard solar model |J, confirming that the electron neutri- 
nos are oscillating into some combination of fi and r neutrinos ||. The signal favors the 
large-mixing-angle (LMA) solution to the solar-neutrino problem, corresponding to mix- 
ing parameters Sm^ aj . ~ 10 _5 eV 2 and sin 2 2# so i ar ~ 0.8. The LMA will be tested by the 
reactor experiment KamLAND 

If all these results are cast in the form of 3-neutrino mixing, then there is still an 
unknown mixing angle, which is currently limited by reactor neutrino data to be 
sin 2 2# roac tor~0.1 [10, 11]. The LSND experiment [12] complicates the picture, requiring a 
singlet neutrino or CPT violation in the neutrino sector fi~3T| . The oscillation parameters 
inferred from the solar and atmospheric results imply that neutrinos will change flavor 
in the outer envelope of a supernova. The details of the transformation depend on a) 
the unknown mixing angle ^reactor and b) whether there are more than three species of 
neutrino, as implied by the combination of atmospheric, solar and LSND results. 

Several existing facilities can detect supernova neutrinos. If a supernova exploded 10 
kpc from the earth with a luminosity of 3 x 10 53 ergs and the energy partitioned equally 
among neutrino flavors, SuperK would see about 8300 events from v e + p — * n + e~ in 
its water detector [|14|]. KamLAND would see about 330 events from the same reaction 
in its scintillator, and SNO would see about 360 such events in its light water component 
KB). KamLAND may be able to measure the spectrum of the high-temperature neutrinos 



through neutral-current neutrino-proton elastic scattering [16]. SNO would also be able 
to see 80 events from v e charged-current break-up of the deuteron. This number would 
increase in the presence of oscillations. An additional 500 events in SNO would come from 
the neutral-current break-up of the deuteron ]fT^1 , which has a low threshold. An analysis 
of the supernova signal in a water-Cherenkov detector and in a heavy-water detector has 
been conducted by in ref. [18]. An analysis of neutrino mass limits from supernova time 
of flight is done in ref. [|H| . 

Here we consider a detector based on, e.g. lead perchlorate (Pb(C104)2), first proposed 
by Elliott ||. Lead had been considered previously as a supernova detector by several 
groups, e.g. OMNIS fl]] and LAND Q. Lead has an attractively large neutrino-scattering 
cross section per nucleon compared with other elements, and most of the scattering events 
produce neutrons [20, 21, 22]. Lead perchlorate, which would be sensitive mainly to the 
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higher-energy neutrinos, has the appealing ability to measure energy deposited by the 
electrons produced in charged-current reactions in coincidence with zero, one, two, or 
more neutrons, as well as to measure the number of neutrons emitted (in isolation) in 
neutral current reactions. Per kt of lead perchlorate, for the supernova described above, 
there would be about 378 electrons produced with one neutron and about 234 with two 
neutrons, if the oscillations are induced only by the large solar mixing angle and the high- 
energy neutrinos have a Fermi-Dirac distribution with temperature T = 8.0 MeV, and 
effective chemical potential r\ = (see eqn. (|j)). In the neutral-current channel, there 
would be about 105 one-neutron events and about 72 two-neutron events. These numbers 
come from the calculations described below. 

In what follows, we investigate the supernova-neutrino signal in lead perchlorate in 
more detail. We discuss how the electrons produced by the charged-current scattering on 
lead can be used to obtain spectral information on the high-temperature neutrinos. We 
then examine the possibility that by using this information together with a comparison 
of the numbers of charged- and neutral-current events, one could distinguish among the 
oscillation scenarios discussed below. 

In section [2] we discuss the particulars of the neutrino oscillations and what they 
would mean for the spectrum of electron neutrinos coming from a supernova. In section 
H we describe our calculations of the cross sections for neutrino-induced spallation of 
neutrons from lead via the charged and neutral currents. We compare our results to 
previous calculations and discuss uncertainties. In section || we present the results of our 
calculations and show how they may be used to learn about the spectra and oscillation of 
supernova neutrinos. Section || is a conclusion. 



2 Oscillation Scenarios 



Although great strides have been made recently in understanding neutrino oscillations, 
there are still unknown entries in the neutrino mixing matrix. The mixing matrix can be 
written as 
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where c 



COS 



9ij and s^j 



smi 



are the rotation angles, the subscripts 1,2,3 denote 
the mass eigenstates of the neutrinos, and 5 is the CP-violating phase. In combining the 
solar and atmospheric data into a single 3 by 3 mixing matrix, one typically associates 
5m%i = m 2 . — m 2 « 5m 2 ^ ar . ~ 10~ 5 eV 2 . The mass-squared difference associated with 

w 3 x 10~ 3 eV 2 . The relation 5m 2 = makes 
eV^. What is not clear is the sign of $777.23, si nce the vac- 
uum oscillation of atmospheric neutrinos is independent of the sign. As a consequence the 
sign of is unknown as well. This has important implications for supernova neutrino 
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oscillations, since, depending on the last sign, either electron neutrinos or electron antineu- 
trinos will undergo a resonance associated with 5m 2 1 . The MSW effect, unlike vacuum 
oscillations, depends on the sign of the mass difference. The sign of Sm^i is determined 
by the requirement that matter-enhanced oscillations in the sun cause a resonance in the 
v e rather than the v e channel. 

By finding the vacuum survival probabilities of electron and muon neutrinos, one can 
further associate sin 2 2# a t mos = sin 2 2#23 cos 4 #13 ~ sin 2 2#23 and translate the limit on the 
electron neutrino mixing angle from reactor experiments to sin 2 2^13 < 0.1. Finally, one 
associates sin 2 26> so i ar with sin 2 1Q\i- 

What are the consequences of the above analysis for supernova neutrinos? In the 3 by 



3 mixing scheme there are two densities at which resonances can occur [23]. In some cases 



shock propagation [24] and nonaxisymmetric explosions [25] will cause the neutrinos to 
encounter their resonance density more than once and thus cause time-dependent features 
in the profile. Here we consider the problem in the first couple of seconds after the 
explosion, before the shock has had a chance to propagate to the point where the neutrinos 
may undergo a particular resonance twice. 

At these relatively early times, there are two main resonances through which the 
electron neutrinos will propagate. The densities at which the resonances occur are set 
by the 5m 2, s and the energies of the neutrinos. For a given density, the energy that a 
neutrino must have to undergo a resonant transformation is given by 

„ , . 5m 2 cos 29 v , , 

Eres(r) EE ± y ^ \ (2) 

where the potential is V(r) = 2^/2G F [N e -(r) - N e +(r)}. 

The 5m 2 3 resonance: The sign of <5m 2 olar is fixed by the solution to the solar neutrino 
problem. However, as noted above, the sign of Sm 2 3 is still unknown. In a "normal 
hierarchy" (5m 2 3 > 0), electron neutrinos, as opposed to antineutrinos, would undergo a 
resonance at a density of about 10 3 gcm~ 3 . (This is for a 40 MeV neutrino and a mass 
scale 5ml 3 m 3 x 10~ 3 eV 2 .) If the resonance were adiabatic, which would be the case if 
sin 2 2#i3 3> 10 -3 (see fig. |]), then the original electron neutrinos would become mostly 
which would be measured in the detector as part Un and part v T . The original v e would 
therefore be seen primarily in the neutral current events. The electron neutrinos would 
almost all acquire the high-temperature spectrum of the original v„ or v T . The v 3 would 
not mix with the other species again in the supernova, so the electron neutrinos would 
retain the "hot" spectrum as they arrive at the earth. 

The 5m 2 12 resonance : If the electron neutrinos didn't transform at the density associ- 
ated with Sm 2 3 , which would imply sin 2 2^13 <C 10~ 3 , they would have a second oppor- 
tunity at a lower density, ~ lgcm -3 , which is associated with <5m 2 olar ~ 10 _5 eV 2 . We 
assume that this resonance is adiabatic, as implied by the presupernova model profiles of 
(|26|| . Then the electron neutrinos as they would be measured as they arrive at the earth 
would have sin 2 6* so i ar of the original v e spectrum together with cos 2 6> so i ar of the original 
Vu,t spectrum. The number of charged-current events in the detector would therefore be 
sin 2 #soiar times the number of events expected with only the lower-temperature spectrum 
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Figure 1: Survival probability as a function of sin 2 2#i3 for an electron neutrino with energy 
10 MeV, 40 MeV and 100 MeV (left to right), as it goes through the resonance associated with 
Smf 3 w 3 x 10 _3 eV 2 . For values below sin 2 6*13 ~ 10~ 3 the neutrinos above 40 MeV pass by 
unaffected, while for values much larger than this, the neutrinos convert. The neutrinos that do 
not convert in this resonance will partially convert at the later resonance associated with Srn^i- 

and cos 2 # S oiar times the number of events expected with the higher-temperature spec- 
trum. The spectrum could be further modified by the earth-matter effect if the neutrinos 
travelled through the earth before reaching the detector. 

We do not consider scenarios intermediate between those just described, i.e. a partially 
adiabatic first resonance, with partial transformation. In other words, as fig. [I] shows, we 
are considering situations where sin 2 29±3 ^> 10~ 3 or sin 2 2#i3 <C 10~ 3 . 

Inverted hierarchy: So far we have assumed a normal hierarchy where the third mass 
eigenstate is the heaviest. An "inverted" hierarchy occurs when this state is the lightest. 
This would cause a resonance to occur in the electron aniineutrino channel at the higher 
density resonance associated with <5m 2 3 . Since a lead-based detector sees very few u e 's, the 
inverted hierarchy would produce nearly the same signal from lead in OMNIS or LAND 
as the normal hierarchy with very small #13. 

We limit our discussion to the situation where there is relatively little earth matter 
effect and are no sterile neutrinos. 
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3 Calculations 



In this work we assume a Fermi-Dirac distribution for the high-energy neutrino component. 
The number flux of neutrinos is given by 

where 

1 E 2 

f{El/) = T^F 2 ( V ) exp(E u /T u - V ) + l (4) 

and (E v ) = T v F^(vi) / '^2(77). The functions F2 and F3 are the Fermi integrals 

Jo exp(x - 77) + 1 

We will assume equipartition of energy between the six neutrinos and antineutrinos, and 
vary the values of T and 77. Because the total energy in each species is fixed, if the average 
neutrino energy rises, the number flux of neutrinos decreases. 

To calculate the event rate from supernova neutrinos we also need to know the cross 
sections of the charged- and neutral-current reactions of interest : 

Bi + n + e" (6) 

Bi + 2n + e" (7) 

Pb + n (8) 

Pb + 2n (9) 

Of course, only half of natural lead is 208 Pb, but the cross sections from the other isotopes 
should not be very different because the thresholds are the same to within a couple of 
MeV and the strength distributions depend largely on collective nuclear motion that is 
not much affected by the addition or removal of a few nucleons. The reaction rates in the 
other isotopes can be calculated through slight modifications to the methods described 
below. 

Even in 208 Pb, these cross sections have not yet been measured. We therefore must 
rely on theoretical predictions. The cross section for a neutrino with energy E v to create 
an electron with energy E e is 

Q2 rl 

a(E u , E e ) = — cos 2 6cPeE e 5{E v - E e - E fi ) J ^ d(cos 9)M f3 , (10) 

where G cos 9c is the effective weak coupling constant, 9 is the angle between the di- 
rections of the incident neutrino and the outgoing electron, p e is the outgoing electron 
momentum, and E^ is the excitation energy of the 208 Bi with respect to the ground state 
of 208 Pb. The quantity Mp is a combination of matrix elements operators of the form 
jl(qr)Yi(tt)T + , ji(qr)Yi(Q)ar + , ji(qr)Yi(tt)VT + , and ji(qr)Y l (Q)aVT + (these include the 
usual Fermi and Gamow- Teller operators in the long-wavelength approximation) where r + 
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changes a neutron into a proton and q = \q\ is the momentum transferred to the nucleus, 
which ends up as 208 Bi. 

Our calculations are based on the coordinate-space Skyrme-Hartree-Fock method in 
a 20 fm box, followed by the A- and .B-matrix version of RPA with the same Skyrme 
interaction in the basis of Hartree-Fock states. We actually do the calculations twice, 
using two different sets of expressions ([^7], ^]) for Mp as well as two different Hartree- 
Fock and RPA codes. For T = 8, r] = neutrinos, the cross sections from the two versions 
differ only by 8% when one neutron is emitted and 4% when two are emitted. We include 
all multipoles up to J = 4. (The contributions of higher multipoles are smaller at these 
neutrino energies.) We truncate the Hartree-Fock spectrum at a value high enough to 
contain most of the strength; increasing it by a few MeV produces almost no change. In 
particular, the configuration space used is large enough to satisfy both the non-energy- 
weighted and the energy- weighted sum rules. We reduce the value of the axial- vector 
coupling constant gA from 1.26 to 1 in the 1 + channel to account for spin-quenching (this 
may still be a slight overestimate for lead; see ref. ^0|). Finally, we treat the Coulomb 
interaction of the outgoing electron with the final nucleus as in ref. (2^] . (We comment on 
the procedure for the Coulomb corrections below.) 

An amount of energy Efi = E v — E e is transferred to the nucleus, so that a number 
of nuclear states can be excited. Much of the cross section is in giant resonances, e.g. 
the Gamow- Teller (GT) resonance, 15.6 MeV above the ground state of 208 Bi, and the 
dipole and spin-dipole resonances, which make up a broad peak centered about 25 MeV 
above the ground state. It's important for calculations to reproduce these resonances; an 
error of a few MeV can make about 20% difference in the flux-averaged cross sections 
[|22[| . In our RPA calculations we used two Skyrme forces, the venerable SIII [^9| and 
the force SkO' |3(| with the "time-odd" terms in the energy functional explicitly adjusted 
to reproduce GT strength distributions as described in ref. [31|; we refer to the latter as 
SkO+. The interactions both correctly locate the isobaric analog state and GT centroid 
to within an MeV, put the dipole and spin-dipole strength in the right region, and give 
similar results more generally. Fig. |2| shows the total cross sections from the two as a 
function of energy. They are nearly identical except at high energies that contribute little 
to the supernova-neutrino flux. 

For neutral current cross sections, eq. (|To|) still holds provided the outgoing electron 
is replaced by a neutrino or antineutrino (the quantity Mp contains different coefficients, 
and isn't corrected for Coulomb effects). We again quench all matrix elements in the 1 + 
channel by (1/1. 26) 2 , and assume that the nucleon carries no isoscalar (i.e. strange) axial 
current. Our calculation here follows the same steps as in the charged-current channel 
except that we use the like-particle RPA instead of the charge-changing RPA. 

We want to calculate the individual cross sections for the emission of 0, 1, or 2 neutrons. 
Following ref. p0[| , we do so by assuming that two or more neutrons are emitted if the 
excitation energy of the final nucleus is more than 2.2 MeV above two-neutron-emission 
threshold, and one neutron if it's less than that but greater than the 100 keV above one- 
neutron-emission threshold. For high-energy supernova neutrinos, the 1 + excitations in 
charged-current reactions are usually associated with the emission of one neutron and the 
dipole and spin-dipole states with the emission of two. For example, with a neutrino flux 
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Figure 2: Total charged-current cross sections, in fm 2 as a function of neutrino energy, for the 
Skyrme interactions SIII (solid line) and SkO+ (dashed line). 



characterized by T = 8, 77 = — typical predicted values for the hot spectrum — the 1 + 
states contribute 53% of the one-neutron cross section and only 5% of the two neutron 
cross section. The dipole and spin-dipole states (0~, 1~, and 2~) contribute 62% of the 
two-neutron cross section but only 15% for one-neutron emission. 

Our results are nearly identical with those in ref . [^] and are close to those of ref . |2l| ; 
the total cross sections as a function of energy agree to within 20 — 30% in the charged- 
current channel and 15% in the neutral-current channel. They are considerably smaller, 



however than those of ref. [20|, as noted in refs. [21, 22]. Although ref. |2(| uses a different 
nuclear model, the reasons for the difference in cross sections are more basic. Ref. p2f 
pointed to the neglect of terms beyond lowest order in the long-wavelength expansions 
of the multipole operators in Mp, and to a different treatment of Coulomb corrections. 
The former increases the cross sections (compared to ours) dramatically, while the latter 
reduces them somewhat, but still leaves them too large. 

We have treated the Coulomb corrections the same way as in ref. [22]: at each electron 
energy we take the smaller of the results with a simple Fermi function and the effective 
momentum approximation p2[. This procedure almost certainly overestimates the cross 



section a little as well, but not badly j33], |22|| . We also tried the procedure of ref. [20]: 
assuming that the outgoing electron was always in an s-wave and correcting its wave 
function according to the tables of ref. 34], which contain the effects of a Fermi function 
for a spherical charge distribution with nonzero radius. This prescription, which typically 
underestimates Coulomb effects at high electron energies, reduces the total cross section 
at T = 8, 7] = by about 20%. We give our results in tabular form (with the first Coulomb 
prescription and the force SIII) in Table |l|. 
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Table 1: Neutrino cross sections in units of 10 -40 cm 2 as a function of energy (MeV) for 
emission of one and two neutrons, and summed over all decay channels, obtained with 
the Skyrme force SIII. We include the charged-current channel for neutrinos, and the 
neutral-current channel for both neutrinos and antineutrinos. 
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How accurate are our calculations? The good agreement between the two quite differ- 



ent Skyrme forces and between our calculations and those of refs. [21, ^] is encouraging. 



It is possible, however, that all such calculations contain similar mistakes. They are all 
done in the RPA, and they all treat Coulomb effects in a similar way. There may be 
systematic nuclear structure uncertainties, related, e.g., to the issue of quenching of for- 
bidden multipoles in charge-exchange reactions, that no data accurately address. It's hard 
to assert, especially without a more detailed comparison to available data from hadronic 
charge-exchange and electron-scattering experiments, that the cross sections are accurate 
to better than 25% for high-temperature supernova-neutrinos. The average recoil energy 
is a much more robust quantity, however; there our two calculations typically differ by 
less than a few percent, and even when we change the treatment of Coulomb effects, the 
difference is typically less than 5 or 10%. We can probably take that range as a reasonable 
estimate of the uncertainty in the average energy. We will make use of this observable in 
the next section to show how the temperature of the supernova neutrinos can be signifi- 
cantly constrained. 



4 Results 

Here we use the calculations described in the previous sections to show: 1) how to use 
charged-current events in which an electron is emitted in coincidence with one neutron or 
with two neutrons to characterize the neutrino spectrum, and 2) how to use the charged- 
and neutral-current signal to determine what sort of oscillation has taken place. 

First, we display in fig. || the numbers of events associated with various parameteriza- 
tions of Fermi-Dirac spectra. Next we consider the electron energies. It takes a neutrino 
of about 18 MeV to excite lead with much probability. The observed electron energies 
are therefore chiefly sensitive to the high-energy tail of the neutrino spectrum. In fig. 
U we show the electron-energy spectra associated with one and two neutrons, obtained 
by folding the cross sections calculated with the Skyrme force SIII (which we use below 
unless we explicitly state otherwise) with the neutrino flux in eqs. (^-||) for T = 8 MeV, 
rj = 0. The spectra observed in the detector will have very nearly these shapes with either 
partial or complete oscillations. The high thresholds associated with one- and (especially) 
two-neutron emission make the signal even less sensitive to the original low-temperature 
component of the spectrum. 

If there were no oscillations at all, and all the electron neutrinos had a temperature 
T = 3 MeV and an effective chemical potential r^ff = 3, there would be only ~ 92 one- 
neutron events and ~ 5 two-neutron events. If oscillations were to occur, these numbers 
would be much larger. If the transformation associated with the <5mf 3 were complete, there 
would be ~ 460 one-neutron events, with an average electron energy of 22.9 MeV, and 
~ 310 two-neutron events, with an average energy of 26.2 MeV. If instead the oscillations 
were due to the resonance associated with <5mf 2 (sin 2 2# so i ar 0.75), the ~ 370 one- 
neutron events would have an energy of 22.1 MeV, hardly any different than before, and 
the two- neutron events would simply be reduced in number by almost exactly 25%, with 
little change in the spectral shape. 
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Figure 3: These contour plots show the number of electrons emitted from one kilotonne of lead 
with one neutron (left panel) and two neutrons (right panel). The contours on the left panel are 
from left to right 500, 450, 400, and 350. The contours on the right panel are from left to right 
450, 350, 250, and 150. The horizontal axis is the temperature of a Fermi-Dirac distribution of 
neutrinos, while the vertical axis shows the effective chemical potential. 



Because the shapes of the electron-energy spectra do not depend significantly on 
whether the oscillations are complete or partial, they can provide a lot of information 
on the original /x/r- neutrino spectrum. Of course, we do not know the cross sections 
perfectly, so extracting the neutrino spectrum from the electron spectra is not trivial. We 
therefore focus on an observable that proved fairly robust against changes in the calcu- 
lation: the average electron energy for one- and two-neutron events. In fig. ||| we show a 
few contours in a plot of the average energy as a function of the temperature and effective 
chemical potential, assuming a Fermi-Dirac distribution for the higher-energy component 
of the neutrinos. The figure shows that the higher the temperature of the neutrinos, the 
higher the average electron energy. The same is true of the chemical potential, particularly 
when it reaches 2 or 3. 

To show that the calculations of average energies are not very sensitive, fig. || presents 
the results with the two different forces, SIII and SKO+. Fig. || shows the same plot 
but with more contours and only for the force SIII. The contour lines in fig. ^ show the 
effects of uncertainties in the calculated average energies. The figure shows the ±5% and 
±10% regions for an electron spectrum with an average energy of 20 MeV emitted with one 
neutron, and for a spectrum with average energy 23 MeV emitted with two neutrons. Even 
with an uncertainty of ±10% (see the last section), the temperature can be determined to 
within about 1 MeV. If the average energies can be accurately measured, therefore, the 
temperature of the neutrino distribution can be determined well. On the other hand, if 
one considers only the average energies, a large range of chemical potentials is allowed (at 
least for r] < 3), almost independent of how well the average energy can be measured. 
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Figure 4: Energy spectra of electrons associated with one-neutron (solid line), and two-neutron 
(long dashed line) events obtained by folding the cross sections (Eq. with the neutrino fluxes 
(Eq. ||) for a neutrino spectrum parameterized by T = 8 MeV and r\ — in one kilotonne (10 6 kg) 
of Pb. 

One might ask how much we gain by considering the one- and two-neutron contours 
simultaneously. We could plot the two sets on the same figure, creating a grid that would 
in theory allow us to determine both T and rj. Unfortunately the contours for the two 
kinds of events are nearly parallel except at large rj, so that the such a plot may not 
provide a stronger constraint on the flux parameters. 

In fig. ^, we plot the width of the distribution, defined as (E^ — (E^) 2 ) 1 / 2 instead of the 
average energy. We use ±5% and ±10% contours for a width of 12 MeV for both electrons 
emitted with one neutron and electrons emitted with two neutrons. This quantity, too, is 
not sensitive to rj unless rj is large. The ±5% and ±10% contours depend, of course, on 
what average energy is measured. Furthermore, a similar analysis, both here and for the 
average energies, may be performed with deviations from a Fermi-Dirac spectrum. 

The extent to which the electron spectra are useful for determining the high-temperature 
component of the spectrum depends very much on the contamination from other neutrino 
sources. As discussed above, in the case of partial transformation, the low-temperature 
component of the spectrum will contribute very little to the average energies of electrons 
associated with two neutrons, but can make almost an MeV difference in the average en- 
ergy of electrons associated with one neutron. Furthermore, if the lead perchlorate were 
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Figure 5: These contour plots show the average energies of electrons emitted from lead with one 
neutron (left panel) and two neutrons (right panel). The contours on the left panel are from left 
to right 24 MeV, 22 MeV, 20 MeV, 18 MeV, and 16 MeV while those on the right panel are from 
left to right 26 MeV, 24 MeV, 22 MeV, and 20 MeV. 

mixed with water, interactions with Chlorine, Oxygen and Hydrogen would also occur. 
Oxygen and Chlorine have high thresholds for the emission of neutrons, and they also 
emit protons, so their event rates will be low. But if the water were 20% of the detector 
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Figure 6: These contour plots show the average energies of electrons emitted from lead with one 
neutron (left panel) and two neutrons (right panel). The two contours on the left (right) side of 
each panel are for 23 Mev (19 MeV) average energy for two different forces. 
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Figure 7: Contour plot showing the average energies of electrons emitted from lead with one and 
two neutrons. In the plot on the right hand side, the bands show the ±5% (inner) and ±10% 
(outer) bands for an average energy of 20 MeV of an electron emitted with two neutrons. In the 
plot on the left hand side, the bands show the ±5% (inner) and ±10% (outer) bands for an average 
energy of 23 MeV of an electron emitted with two neutrons. 
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Figure 8: Same as Fig. 0, except that here we show contours of the width of electron energy 
distribution for one (left figure) and two (right figure) neutrons. We show contours of ±5% and 
±10% for a width of 12 MeV. 



there would be approximately 70 events from 

v e + p — ► n + e 



fill 



(assuming a normal hierarchy with sin 2 #12 ~ 0.75) per kt of lead perchlorate. These 



14 



events will be difficult to distinguish from the one-neutron charged-current events on lead. 

In the case of a normal hierarchy and no earth matter effect, the antineutrino signal 
would consist of ~ sin 2 #12 of the original v e spectrum and cos 2 #12 of the original v^^ 
spectrum. If a detector such as SuperKamiokande or KamLAND were operating at the 
time, it would be able to determine the v e spectrum, making it possible to subtract out that 
component. If not, however, one might have to rely on the two-neutron signal. That by 
itself, though, is enough to get a good handle on the temperature of the high-temperature 
neutrinos, as discussed above. 

We turn finally to the problem of distinguishing between two likely neutrino spectra at 
the detector: a completely hot spectrum associated with the normal hierarchy and com- 
plete conversion at the first resonance, and a ~ 75% hot, ~ 25% cold mixture associated 
with either the inverted hierarchy or very small #13. Once the electron energy spectrum 
has been used to learn about the hot component of the electron-neutrino energy spectrum, 
we can then use the ratio of charged-current to neutral-current events to get around the 
problem of not knowing the distance to the supernova. 

Fig. P shows the ratio of the number of charged-current events to the number of 
neutral-current events (with no electrons). Since in the charged-current reactions the 
neutron emission threshold is high, only high energy neutrinos cause events. Furthermore, 
the neutral current signal is unaffected by oscillations. Therefore, the partial-oscillation 
scenario results in a ratio that is about 3/4 ( or cos 2 #12) of what complete oscillations 
yield. Clearly, to exploit this result it is essential to know the reaction cross sections 
precisely, i.e. to within about 10%. No calculations, including ours, are that certain but 
improvements in the model, e.g. through an exact treatment of Coulomb effects, and 
(especially) measurements of neutrino/lead cross sections could get us to that level. 

If experimental work on solar neutrinos continues, then at the time a supernova ex- 
plodes there may well be a detector in place that can see electron antineutrinos through 
the reaction v e +p — > n + e + . The positron spectrum, which will closely track the incoming 
antineutrino distribution, could in principle be used in conjunction with lead to determine 
the difference between the normal and inverted hierarchies. 



5 Conclusions 

In this paper, we have analyzed the properties of a lead-based supernova-neutrino detector. 
The expected energy distribution of electrons emitted with one or two neutrons peaks 
somewhere in the low 20's of MeV, depending on the details of the incoming neutrino 
spectral shape. We expect a few hundred charged- and neutral-current events per kt of 
lead. 

We have used the Random-Phase- Approximation with effective Skyrme forces to cal- 



culate neutrino-lead cross sections. Our results are in agreement with those of refs. [21] 



and [22 1, but lower than the ones in ref. [20|, for reasons we understand. We used two dif- 
ferent Skyrme forces, SIII and a parameterization we call SkO+ based on the force SkO'; 
the average electron energies obtained with these two interactions are within a few per- 
cent of each other. Despite this good agreement, it is difficult to assess quantitatively the 
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Figure 9: Ratio of charged current events to neutral current events in the case of the various 
oscillation scenarios. From left to right: total conversion at the 8m\ s resonance, partial (adiabatic) 
conversion at the 5m\ 2 resonance, and no conversion. This figure is made assuming T=8.0 MeV 



r\ = 3 for the high-temperature flux, T=5.0 MeV, r\ 
i] = 3 for the original v e flux. 



for the original v e flux, and T=3 MeV 



overall systematic uncertainties in our cross sections. We don't know how much forbidden 
strength, which is particularly important in the two-neutron channel, is quenched, or even 
much about how it's distributed. It's possible, however, to improve the calculations, and 
additional data from neutrino scattering and electron scattering would help tremendously. 

The spectra of neutrinos reaching the earth can be modified by oscillations. We have 
argued that in most scenarios, the electron energy spectrum can be used to determine the 
temperature (and perhaps the effective chemical potential) of the hot neutrino spectrum, 
due to neutrinos that were originally emitted as or v T . (We find that lead has little 
sensitivity to the cold neutrinos originally emitted as v e , or, in the charged-current channel, 
to antineutrinos.) With a ±5% theoretical uncertainty in the average energies of the 
electrons, there would be an approximately 0.5 MeV uncertainty in the temperature of 
the hot neutrino spectrum. This doubles if the uncertainty is ±10%. These numbers 
are small enough so that a lead detector would provide crucial information about the 
proto-neutron star. 
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We discussed the possibility of distinguishing between the complete and partial neu- 
trino transformation by using the ratio of charged-current to neutral-current event num- 
bers. Although the idea works in principle, it would be hard to use now because of the 
uncertainty in calculated cross sections. However, it should be possible to distinguish be- 
tween the case of no transformation and either partial or complete transformation. Since 
some degree of transformation is expected, either due to the LMA or #13 this ability will 
provide an important check of our understanding of neutrino physics. 

A measurement of neutrino cross sections on lead is very important. As discussed 
above, we need a reasonable degree of certainty in the average energies to extract infor- 
mation about the incoming neutrino spectrum. That certainty is even more important 
if we want to see whether some or all of the original electron neutrinos have oscillated. 
The ability to distinguish between the two would let us use supernovae to learn some- 
thing about the size of #13 or discover whether nature has chosen the normal or inverted 
hierarchy. 
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